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Abstract 

The theory of the Clausius' virial maximum to explain the Fundamental Plane 
(FP) proposed by Secco (2000, 2001,2005) is based on the existence of a maximum 
in the Clausius' Virial (CV) potential energy of a early type galaxy (ETC) stellar 
component when it is completely embedded inside a dark matter (DM) halo. At the 
first order approximation the theory was developed by modeling the two-components 
with two cored power-law density profiles. An higher level of approximation is now 
taken into account by developing the same theory when the stellar component is 
modeled by a King-model with a cut-off. Even if the DM halo density remains a 
cored power law the inner component is now more realistic for the ETGs. The new 
formulation allows us to understand more deeply what is the dynamical reason of the 
FP tilt and in general how the CV theory may really be the engine to produce the FP 
main features. The degeneracy of FP in respect to the initial density perturbation 
spectrum may be now full understood in a CDM cosmological scenario. A possible 
way to compare the FPs predicted by the theory with those obtained by observations 
is also exemplified. 
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1 On the tilt 

It is well known that galaxies of different morphological types cluster around 
the Fundamental Plane (FP) (Dressier et al., 1987; Djorgovski and Davies, 
1987; Faber et al., 1987; Bender et al., 1992; Djorgovski & Santiago, 1993; 
Renzini & Ciotti, 1993; Ciotti et al. 1996; J0rgensen, 1999; see, e.g., the review 
of D'Onofrio et al., 2006, and the references therein) in the three dimensional 
space of: Ve, effective radius; Je, the mean effective surface brigtness within 
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re] (Toi the central projected velocity dispersion. On the basis of homology + 
virial theorem one would expect that the FP equation has to be: ~ cf^I^ 
where A = 2, B = —1. That results completely in disagreement with the 
observations in different bands. Typical values in 5— band are: A = 1.33±0.05; 
B = —0.83 ± 0.03 (e.g., in D'Onofrio et al., 2006). These unexpected values 
produce in the k coordinate system (Bender et al.,1992) the so called tilt that 
is an increasing of the ratio: dynamical mass M^yn over luminosity L, of this 
kind: 

MdyjL ~ (M,,„)°-2 (1) 

Many attempts have been done in order to understand the FP tilt which is 
also one of the common features either for galaxy FPs or for the FPs of all 
virialized structures which all together define the so called cosmic meta-plane 
(Burstein et al., 1997). The review of D'Onofrio et al. (2006) may help the 
reader to take into account the more recent efforts to solve the hard problem of 
finding an explanation of the trend ([1]) when the K — band is also considered 
and then the population effect has to be rouled out. Actually it is possible 
to explain the trend observed in the i?— band as a metallicity sequence of an 
old stellar population (Maraston, 1999). However the Mdyn/L values in the 
/T— band are independent of metallicity even if the tilt is observed (Pahre 
et al., 1998). A secondary effect is then needed to explain the iiT— band tilt 
(Gerhard et al., 2001). 

The Clausius' Virial theory (TCV) of FP has the aim to propose a dynamical 
mechanism able to produce the required effect on a huge range of mass scales 
from globular clusters to galaxy clusters. The purpose is to prove that it 
may be possible to change A, B exponents (from the expected values 2, —1) 
without breaking homology + virial equilibrium. It is based on the existence 
of a special virial configuration characterized by a maximum in the Clausius' 
Virial potential energy (CV) which, on galaxy mass scale, refers to a baryonic 
(stellar, B) component when it is completely embedded inside a DM halo {D 
component). At the first order approximation (linear) the two-components 
are modelled with two power-law density profiles and two infinitesimal cores. 
The general strategy is described in many papers (Secco, 2000; Secco, 2001, 
hereafter LSI; Marmo & Secco, 2003; Secco, 2005, hereafter LS5). 

Now we move from a linear approach of TCV to a non-linear one that is to 
an higher level of approximation in which the stellar component is built up 
by a King-model with a cut-off. Even if the DM halo density remains a power 
law the inner component is actually more realistic for the ETGs. The new 
formulation allows us to understand more deeply the physical reasons which 
produce the FP tilt and the role of the main involved quantities, particularly 
that of le- Moreover we may begin the comparison between the expected edge- 
on FPs with those obtained by observations (e.g. that of Djorgovski & Davies 
(1987)) and try to reproduce in K-space the tilt fit-equation of Burstein et 
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al. (1997). Its theoretical derivation may explain why the FP is degenerate 
in respect to the initial density perturbation spectrum in a CDM scenario as 
already underlined by Djorgovski (1992). Some initial examples of theoretical 
FP calibration will be given in the sects. 7, 9, for some special choice of 
theoretical parameters. A more complete discussion is still in progress. 



2 General strategy of TCV 



Briefly summarizing, the general strategy consists to use the two-component 
tensor virial theorem (e.g., Brosche et al., 1983; Caimmi & Secco, 1992) to 
describe the virial configuration of the baryonic component embedded in a 
DM halo at the end of relaxation phase (see, Bindoni & Secco, 2008). It reads: 

2{T,)ij = (K)i,-; (u = B, D- i, j = X, y, z) (2) 

According to the scalar virial for one component, the potential energy tensor, 
which has to enter into the tensor virial equations, is the Clausius' virial 

tensor, (V^).y, built-up of the self potential- energy tensor, {flu)ij, and the tidal 
potential- energy tensor, {Vuv)ij- Then, according to the scalar virial theorem, 
the trace of CV tensor, in the case of stellar component, has to be read: 



VB^nB + Vbd (3) 

r=l dXr 

3 9$, 



(Vbd) = J PsYl^Xr-^ daf^ = J pB{fh ■ fo) dx'h; 



where is the B component density and fs, fo are the force per unit mass 
due to the self and DM gravity, respectively, at the point rg and are 
the respective potentials. 

Conversely, the total potential energy tensor of the B component is: {^B)ij + 
(W^BD)ijj where the interaction energy tensor is: {WBD)ij = ~\ J PBi^D)ij dxB', 
and the potential tensor due to the DM (e.g., Chandrasekhar, 1969) is: {^D)ij = 

GJpn{x')^-^^^^^^^P^dxn. 

\x—x\ 

To be noted that in general: (Vbd)^ (W^BD)jj , the difference gives the resid- 
ual energy tensor (Caimmi & Secco, 1992). 

We will describe a re-formulation of TCV in the case in which the two- 
component model is built up of: a bright B stellar component with a King 
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(1962) truncated density profile completely embedded in a DM frozen halo, 
D, with a cored power-law mass density distribution. 



3 Why introducing King's model 



3. 1 End of relaxation phase 



The violent relaxation mechanism leads to an equipartition of energy per unit 
mass and not per particles (see, e.g., the review of Bindoni & Secco, 2008, 
and references therein). If a is the velocity dispersion, assumed to be the same 
for every star mass, integration of the distribution function, f{E), over the 
velocities (Binney & Tremaine, 1987, Chapter 4; Combes, 1995, Chapter 4), 
yields the density: 

p(r) = pie-^(^)/-' (4) 

where the total energy per unit mass is: E = {l/2)v'^ + U; [v and U are velocity 
and potential energy per unit mass, respectively). On the other hand, Poisson 
equation: 



-,^{r'^)=AnG ff{E)dv (5) 
ar ar J 

becomes by means of Eq.(j4]): 



with the solution: 



In turn, Eq. (jll) gives: 



^2 



- 2^ 



2H^^) = U{r)/a' (8) 

when a core radius rc = 3a{ATTGpo)~^^'^ is introduced in order to avoid an 
infinite value of the central density po- Tc corresponds to the radius at which 
the projected density of the isothermal sphere falls to roughly half of its central 
value. Eq.([8]) gives us the asymptotic behavior as soon as r is greater of about 

U{r) ^2(r^\n{r/rc) (9) 
which means again from Eq.(jll), an isothermal behavior, p(r) oc r"^ as r — > oo. 
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3.2 Problems with isothermals 



The isothermal energy distribution function extends spatially to infinity with 
infinite mass and so does not be suitable to represent a real elliptical galaxy. 

Since 1965 Ogorodnikov has highlighted that: in order to find the most prob- 
able phase distribution function for a stellar system in a stationary state, the 
phase volume has to be truncated in both coordinate and velocity space. While 
in the velocity space the truncation arises spontaneously due to the existence 
of escape velocity, the introduction of a cut-off in the coordinate space ap- 
pears, on one side, necessary in order to obtain a finite mass M and radius R, 
but, on the other, very problematic. 

A similar difficulty also appears on the thermodynamical side, for which an 
extensive literature exists (from: Lynden-Bell & Wood, 1968; Horowitz & Katz, 
1978; White & Narayan, 1987, until, e.g., Bertin & Trenti, 2003, and references 
therein). By using the standard Boltzmann-Gibbs entropy: 



defined by the distribution function, f{x,v) (hereafter DF), in the fi phase- 
space, and looking for what maximizes the entropy of the same stellar system, 
the conclusion is: the DF which plays this role in (fTOj) is that of the isothermal 
sphere. But, the maximization of S, subject to fixed mass M and energy E, 
leads again to a DF that is incompatible with finite M and E (see, e.g., Binney 
& Tremaine, 1987, Chapter 4; Merritt 1999, Lima Neto et al. 1999, Marquez 
et al. 2001, and references therein). 

Our limited contribution to the wide discussion existing in the literature will 
be to underline as in a stellar component, embedded in a second dark matter 
subsystem (e.g., Ciotti, 1999, and references therein), a truncation is sponta- 
neously introduced in coordinate space, due to the presence of a scale length 
induced from the dark halo, as long as virial equilibrium holds. That is the 
tidal radius which has been discovered in the TCV dynamical theory (LSI, 
LS5) we will revisite in the next paragraphs^ 



Even if some considerations which follow are more general and may also be ex- 
tended to spirals, we will limit our considerations to the collisionless stellar systems, 
as the ellipticals are considered. 




(10) 
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3.3 King's models with cut-off 



We will assume for the B-component the empirical surface light density law 
proposed by King (1962) as: 

^^^^ ^ ^4 [1 + (R/Rcn/' " [1 + {Rt/Rcyv/^ ] ^^^^ 

where Rt is the value of R at which I reaches zero. The law has the advantage 
to take into account the existence of a cut-off in the surface density distribution 
as one expects for globular clusters (GC) (the profile ( fTTI) is born for them) 
but also for ellipticals too. Indeed we will need of a similar truncation radius 
because in the TCV theory the meaning of the radius at the special Clausius 
virial maximum configuration has an analogous role (LS5) of that discovered 
by von Hoerner (1958) for GC. The last one is due to the tidal effect of 
Galaxy, the former one to a similar tidal effect but, this time, due to the 
dynamical effect of DM halo distribution on the baryonic component. The 
value Rc corresponds to the core-radius and the value is linked to the 
central surface light density Jq by: 



1 



[1 + {Rt/Rcn/' 



:i2) 



4 King's model in phase-space 



From the previous considerations, an acceptable distribution function in the 
phase-space has to have a cut-off at the energy Eo as King introduced in his 
model (1966): 



fK{E) = OforE>E, (13) 
fKiE) = (27ra2)-3/2p^(e(^°-^)/'^' - 1) for E < E, (14) 

The spatial density can be obtained after integration over the velocities in the 
following way: 




erf being the error function. 
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Following King (1962), the spatial density is given by: 



PB [.z) 
where 



3 r- 



COS z 



zV- 



(16) 



1 + {r/r^f 
.1 + {n/rcY 



.1 + {rt/rcY 



(17) 



with Tc {=Rc) and rt {=Rt) are the core and the cut-off spatial radius, respec- 
tively. The mass inside z is given by: 



M{z) = A'KrlkMZo 







1 r 


f 


[--ll 




J Zo 


. o 


z - 



-l~_(l_^2^1/2 



COS z 



dz 



(18) 



The projected density is: 



S(i?) = k 



M 



(19) 



which is linked to the spatial density by the Abel integral equation (Binney 
& Tremaine, 1987, Chapt.4): 



p(r) 



1 r* rfS dR 



71 Jr dRy/W 



(20) 



S(i?) becomes the classical surface luminosity density I{R) given by King 
(1962) (Eq lTTl) as soon as kM translates into ki (see, FigH]). It should be 
noted that in our models kM ^ ki being one of our main assumptions. Indeed 
in the TCV the galaxy homology family is intrinsically characterized by a ratio 
mass/luminosity for the i? component different from a constant. 

Integration of S(i?) with respect of 27rRdR gives the total projected mass 
within the projected distance R of the center which becomes the luminosity 
function, L{X), with the substitution of S(-R) — > I{R) and kM — ^ ki- 



L{X)=TxrlkLFL{X) 



(21) 



where: 
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ClO 

o 




og(R,/rJ = infty 



log(R/r^ 



Fig. 1. The standard curves of King's model of I{R) normalized to the central value 
as function of the different parameters Rt/vc- The limit case when the cut-off goes 
to oo is also shown (King, 1962). 

Fl{X)= ln(l + X)- + ^ 



1 



X 



X 



R 



1/2 



X, 



l + Xt 



22 



(23) 



According to Eq.(!2Tl). as Xt » 1 the limit of L{Xt) goes approximately to: 



L (Xt) - TirlkL In 



20r2 



(24) 



It should be noted that, according to Eqs. (fT8l [2T!) . the following equation has 
to hold: 
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M(l)^ 



L{Xt 



Mb/L = kM/kL 



(25) 



The trends of normaUzed L{R/rc) and M{r/rc) are shown in the Fig.( 
the case of rt/vc = 20, by assuming kM = ki- 



rn 



L(R/rc) 
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M(r/rc) 




II r/r^ 

5 10 15 20 R/j- 

Fig. 2. Trends of ^^^t^"^ and ^^''A"'* in the case of, rt/vc = 20, by assuming k^j = k^. 
The upper curve is corresponding to the integration of surface luminosity density 
and the bottom one to the spatial mass density integration. At Rt = rt = 20rc the 
Eq.(l25]) holds. 



5 Inserting King's model into TCV 



The aim is to translate the hnear formulation of TCV into a non-linear one 
modeling the bright (or baryonic) inner component, B, of mass Mb with a 
King's model and the dark matter halo, D, of mass Md with a cored power- 
law. 

The first step is to build up the ClausiusVirial tensor trace of Eq.(I3]) which 
for similar strata spheroids is given by: 

Vb = -i^nB -F - vy (26) 

a-B clb 

where z/^s and, uy, are the self mass distribution and the interaction coeffi- 
cient, respectively. F is a form factor which is equal 2 in the spherical case we 
choose for the sake of simplicity. 
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We will follow the general method proposed by Caimmi (1993). Then we define 
(Roberts, 1962; Caimmi & Marmo, 2003, and references therein): 



Fui^u) = 2 /' mu)^ud^u ;u = B,D 



where fuiCu) are the dimensionless density profiles with: 

Cu, o,u, ^ou being the concentration, virial and scale radius, respectively. For 
the King's inner component qb is the cutoff radius rt, Tob — fci Cb — Tt/vc- 

The King's profile, normalized to the scale radius density value, omiting for 
sake of simplicity some obvious indices, is: 



/b(0 = Y 



H 



1/2 



cos 



1 + C2 



(27) 



2\ 1/2' 





\c' - {cm 






1 + C2 


l/2j 



l + C^ 



• COS 



-1 



1 + C2 



(28) 



The power-law profile for DM, normalized again to Po(C — ^/Cd), becomes: 

2 



Now we can define all the coefficients we need for the Clausius trace: 



(29) 



9 1 

yv{x) = --7 — 7 — — — mw^''^*^(x) aB/ao 

o [yBjMXyDjM 

3 

{^u)m = - Fu{iu)diu ■,u = B,D ; 
z Jo 



and: 



9 

The mass ratio Dark/Bright is: 

Md 



m — 



(30) 

(31) 
(32) 

(33) 
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and the interaction term inside the tidal tensor trace, Vbd, (via the coefficient 
z/y of Eg. (1301) ) due to dynamical effect of DM on the baryonic one, is given 
by: 

^(-*)(x) = r FB{iB)'^iDdiD ; iB = in/x (34) 

where {vu)m is an additional profile factor which gives the mass of the two 
components: 

47r 

Mu = {vu)mMou] Mou = yPo^a' (35) 
The main functions cited here are explicitely given in Appendix. 



6 Main features of TCV revisited 

Before to draw the main lines of TCV let us to depict the cosmological envi- 
ronment inside which the approach, based on tensor virial theorem extended 
to two components, tries to explain the galaxy FP- tilt. As we have showed in 
LSI, it is impossible to give account of the galaxy scaling relationships without 
considering the cosmic scenario even if the FP as a whole appears to contain 
a degeneracy in respect to the initial density perturbation spectrum. That has 
been pointed out ffist by Djorgovski (1992). We try to explain why in the sect. 
9. 

6. 1 Cosmological framework 

Our framework is a hierarchical CDM scenario (see, e.g.. Coles & Lucchin, 
1995). We assume that the spherically-averaged properties of a galaxy assem- 
bled dark halo of mass M£, formed by hierarchical clustering, may be deduced 
from the linear theory where the mass variance, cta/^, in a random Gaussian 
field of an Einstein-de Sitter model (Silk, 1999, Chap. 3), evolves from the 
recombination time tree forwards, as : 

aM,(t)~M^("+^)/^(t/t_)2/=^ (36) 

where n = Uree is the effective index. This is a self-similar toy model in which, 
if the linear regime ends at maximum expansion time, tmax (i-e., the free-fall 
time Tff) of the spherical top-hat filtered mini- universe of comoving radius 
aD, mass Md and density po, it holds: Pd]/"^ ~ tjj ~ (MD/a^^)"-^/^, where 
(^D = Ci£)o and Pd = pco at tmax- From Eq. fl5^ it follows that: 

^//~^Af^(W)-'/'~M|" (37) 



11 



where the local slope: 



a = arec = TT^TT = i^rec + 3)/6 (38) 

amMD 



and then, in turn: 



ano ~ M^/^', l/y = (5 + n,,,)/6 = (39) 



The "formation" of halo, i.e., its virialization occurs at Sp/p = 1.67 when the 
time from tree is about tp = 2Tff and its radius becomes: 

an ^ aDo/2 ~ M^/^' (40) 

That is because the collisionless halo system relaxes under the violent re- 
laxation mechanism by conserving its total energy. This is also true even if 
baryons and dark matter particles relaxe together at least under some condi- 
tions (see, e.g., LSI, parag.8). Moreover, at tp the mean density pp is approx- 
imately ISOpu where pu = 1 / (GirGt'j^) . Due to the scaling of tp ^ M^^'^^^^^ 
proved before, at virial equilibrium it holds: pp ~ soon as the 

halo density profile is of Zhao (1996) family, from Eq. fl5^ the mean density 
is given by {i'd)mPoD and then if the concentration Cd, on which {i'd)m de- 
pends, correlates loosely with the mass (Dolag et al., 2004), the central density 
has to scale as: 

Pod - M^("+^)/^ (41) 

Then low mass halos are significantly denser than more massive systems. That 
reflects the higher collapse redshift of small halos (Navarro et al., 1997). Indeed 
the relationship between mass and formation redshift, zp, defined as the z at 
which an object of present mass Mp) has on average acquired half its mass is 
(Lacey & Cole, 1993): 



Zp 



^2(n+3)/3 _ 1)1/2(Mb/M„,)-("+3)/6 (42) 



Mni = 4-lO^'^(l-|-2;)~®/^^^"''M0, is the mass scale over which galaxy count fluc- 
tuations have unit variance, corresponding to the comoving, Rni = 8h~^Mpc. 
These scaling relations are valid provided that the effective spectral index is 
in the range —3 < n < 1. The effective index on typical galaxy scale for a 
scale-invariant initial spectrum is indeed approximately —2 (Gunn, 1987, Silk, 
1999, Chap. 3) and that corresponds to 7' ~ 2. 



6.1.1 Adiabatic contraction factor 

In the primeval version of TCV (LSI) we asssumed that the baryonic com- 
ponent when it reaches virial equilibrium is completely done of a collisionless 
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star fluid. That may be considered as an extreme assumption we have now 
reahstically to soften by a parametrization. We have to assume that when the 
main relaxation process ends the transformation of gas into stars does not 
be completed. To take into account this physical condition we refer to the 
fluid approximation used by Klar & Miicket (2008) (hereafter KM8) in order 
to follow the dynamical coupling of both components the DM one and that 
of baryonic gas. In this approach the gas is allowed to undergo cooling pro- 
cesses and the consequence of its dissipation on DM is an adiabatic contraction 
which, in turn, has a back-reaction onto the gas dynamics. The assumption is 
made that the DM is dynamically reacting fast enough on any change of the 
overall gravitational potential via processes of violent relaxation kind. 

Two cases are taken into account in KM8: i) the initial DM profile follows a 
NFW (1996) profile, ii) the DM distribution is politropic. Even if the approach 
appears very crude for our context (e.g., the whole baryon mass is in gas with 
a ratio of baryons over DM equal to 0.25) we may conclude from both cases 
analized that the effect may be parametrized by introducing in TCV a mean 
linear contraction factor c = fs/f^) < 1/10. 



6.2 Interaction term 



As usually done, the Clausius Virial trace normalized to GM'^F/a^ is given 
by: 



1/. = - ^ (43) 

X X 

Then we have to perform the interaction coefficient uy (Eqs. [301 15^ for dif- 
ferent King's and DM concentrations in the d-range -i- 1. The results are 
collected in the Tables of Appendix together with the explicit formula of Vb- 

In Fig. ([3]) its typical trend is shown in the case: Cb = 10, Cd = 10, d = 
0.5, m = 1 -^20. 

It is to be remarked that also in this non-linear approximation, at good extent, 
the ratio: 

v'y = ^\ ' ~ const. 44 

exactly as in the linear formulation of TCV. That simplifies enormously the 
translation of the linear approximation into the non-linear one. As the first 
consequence is that the Clausius Virial maximum (CVM) appears again at: 
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Vb 



X 




Fig. 3. Normalized Clausius Virial Vb vs. the ratio of virial radii of the two compo- 
nents: X (= aB/cLo)-, for Cb = Co = 10 and m in the range 1 20. 

Moreover M^, which is the fraction of D matter exerting dynamical effect on 
B according to Newton's first theorem, becomes at a good extent, again as in 
the hnear case: 

~ (46) 



The same mass fraction normahzed to Mb becomes: 

m = mx^-'^ = Md/Mb 



(47) 



Owing to Eq. (H5|) , when we consider the special configuration at the maximum, 
the normalized mass fraction is: 



u'y 2-d 



(4J 



which is independent of the mass ratio m. 



Then total mass inside the B-structure which exerts a dynamical effect on B, 
at CVM becomes: 



M, 



dyn 



Mb+Md= MBil+TJlM] 



(49) 



6.3 Energy equipartition 



The presence of Clausius' virial maximum means the virial energy equipartion 
at Xm- It means: 
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- Vbd 



(50) 



Using the definition of masses (Eqj35l) by means of tlieir (z/„)m coefficients 
(1321) . energy equipartition translates into tlie link between tfie two central 
densitieafj as follows: 



PoB ^ 7 — Xm = CXm (51) 

where we have taken into account also the adiabatic contraction (see, KM8, 
subsect.6.1.1) by introduction of the parameter c. 

From the physical point of view Eq. flHT]) means a strict link between the two 
gravitational potential wells of baryons and of DM at the special virial con- 
figuration corresponding to CVM. 



6.4 Light vs. DM halo 



The King's model relationships allow us to link easily the DM potential well 
with the light quantities of the baryonic component in the following way. The 
central mass density of the halo, which defines how deep is the corresponding 
potential well, is linked to Pos (Eq JHT]) . In turn the ratio between the two 
Eq.llISlDSD gives: 



PoB 



S (0) Zo 



Zn) TTTr 



— COS Zo 

Zn 



ZaY 



(52) 



From the other hand Eq.([l9]) reads: S(0) • |^ = /(O). Then by Eqs. ([521 EI]) we 
obtain how the central surface brightness in fiux, /(O), links to DM potential 
well: 



m = F(..)vrr.(l - z^r^^^Pon^P^ 

J^nB {^b)m 

Fizo) = - ^ 



(53) 



— COS~^ Zo 
Zo ° 



2^l/2 



According to Eqs. (l21ll22l) . it is by definition : 



^ Due to the adopted formalism of sect. 5, central density means scale radius density 
value. 



15 



L (Xe) = \7irlkLFL (Xi) (54) 

the solutin of which gives Xg, the square of effective radius normahzed to Tc- 
On the other hand the following relationship for total luminosity holds: 



Uot = 2nr% = nr^hFL (XJ (55) 
and then: 



h = l^k^FL (X,) (56) 



The ratio between Eq. fl56|) and: 

I{0) = kL[l-Zof (57) 
immediately yields: 



I{0) = 2hX,j-^^[l-z,f (58 

Inserting it into Eq.( l53|) we obtain how the quantity of light given by, le 
depends on the DM potential well: 



'^'v {^d)m p ^ 1 ki Fl (X f 
^nB {^b) m ^i-i 2Xe 



Fiz^)nr,^'^Pon^^'-^ (59) 



That is one of the main relationships the TCV yields in order to understand 
the physical tilt-mechanism. We will come back later (sect. 8). 



7 Theoretical FPs 



To find the theoretical FPs in the present non-linear theory approximation (a 
B King model embedded into a D power law halo) becomes easy due to some 
reductions of this approach to the linear one. Indeed as soon as the condition 
dHj) holds the whole main linear formalism (LSI) may be recovered. Then 
there are two ways in order to write down the theoretical equation of PP. 
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>-=Ac~iv 5 10 



Fig. 4. Solution of Eq. (j54p in order to find rg/rc, in the case rt/r^. = 10, of B King's 
component. 

7.1 Main way 



Briefly speaking, the physical reason for the existence of the FP lies into 
the existence of a maximum in the Clausius virial energy which is able to 
divide it in about two equal amounts: the self-potential energy of the baryonic 
component and the tidal potential energy due to the fraction of DM halo which 
has dynamical effect on it. On this special virial configuration the following 
relation holds: 



1 fl^"'^ 

-MB{(y^) = -Vbd = p'vGMbMdF^ (60) 

Extracting om, i-e., the virial dimension of B component at the maximum 
which directly links to re, the FP springs up: 



Factorizing the Eq. (l6T!) as: 



re ~ ao^-^ai-^m-^M^^-^ = (62) 

-T I'Q.B 



it follows immediately: 



Oq — Oo 

3-d , (63) 
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kn, ky are the usual coefficients for kinematic and density galaxy distributions. 
On the contrary of the linear approximation, the present model allows us not 
only to give explicitely the numerical factor in the second relation of fl63p but 
to understand deeply the physical meaning of the previous factorization. 



7.2 The most easy way 



The most easy way to obtain from TCV the theoretical FPs is the following. 
From two-component virial equation (Eq. [2]): 

Tb = -^^B - ^Vbd (64) 

by remembering Eq. (1261) . the definitions of Tb, cr^ and rg in Eqs. ( 1601 1611 1621). 
it follows: _ 

— = kn H kn (65) 

kv re vqb re 

and then Eqs. fl46|48l) . we obtain: 

'2-d 



Mb = cr^r-eCa 



3-d 



(66) 



where C2 = (Gknk^)^^ turns out to be a constant, if we assume that homology 
holds for kinematic and density distributions of elliptical galaxies. Here we 
also assume that C2 comes out from King's model as given by Bender et al. 
(1992, Fig. 5) in the case of isotropic velocity dispersion and with an unchanged 
distribution even if the B King's component is now embedded in a DM halo. 

When Eq.( !66l) is divided by L = cilevl (here ci = 27r), due to the special 

configuration of CVM, then L/Mb ~ Mg^"'*. So the theoretical FP arises in 
the form: 

Te = {c,c,)^cf{L"r^{M%)-^a^lf (67) 

where: 



2 , , 

A = B = 68 

2-d' 2{2-d) ^ ' 

ci = 27r; C2 = 02(^5); C3 = ! (69) 

L° and are luminosity and mass of one elliptical galaxy choosen in order 
to calibrate the plane. 
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1.43 


-0.86 


6.60 


11 


10 


1.43 


-0.86 


6.81 


11 


15 


1.43 


-0.86 


6.94 



Table 1 

Values of parameters which enter into Eq. (j7ip . 
7.3 To test the theoretical FPs 



We compare the theoretical FPs produced by Eq. (l67|) with that obtained by 
Djorgovski & Davies (1987) (hereafter, DD87; see also, Kormendy & Djor- 
govski, 1989) by fitting the observations (in the Lick band): 



log re = 1.39(logao + 0.26(/i)e) - 6.71 
We plot in Figs. ([5]), IQ the edge-on FPs as follows: 



log Te + const 
A 



logao-0.4-(/x). 



(70) 



(71) 



where the values of the parameters are given in Tab.(l). For both figures 
the theoretical FPs from Eq. (!67|) are plotted: i) for d = 0.5 (FigJSj) and 
M%/L° = 15,10,5,1 (from top to down, long-dashed, dot-dashed, dotted 
and short-dashed lines, respectively) and: ii) for d = 0.6 ( FigE]) and ra- 
tios M^/L° = 15, 10,5, 1 (from top to down of, with the same previous line 
types) and the DD87 (solid line) is also shown as comparison. 

To calibrate the theoretical plane fl67|) we use L° and of the elliptical 
galaxy PGC045032 (1300.4 + 2807) of Coma Cluster which has been fitted by 
King's profile (Oemler, 1976) as follows: 



= -16.89; log(rt) = 1.10{Kpc); hg{r^) = 0.10{Kpc) 



(72) 



By conversion of M'^ = — 2.5 log(/(0)r^) + const (in V band) into: /(O) 
3.08 • IO^Lq/IKpc^ and then using Eqs. (I57f58ll5^ of King's model (r^/rc 
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3 

log (pc) 



Fig. 5. Theoretical FPs as outputs of Eq. (j67p for d = 0.5 with different calibra- 
tion M'^/L° = 15, 10, 5, 1 (from top to down, long-dashed, dot-dashed, dotted and 
short-dashed lines, respectively) in re band. The fit-equation of early-type galax- 
ies obtained by Djorgovski & Davies (1987), in the Lick band, is also plotted 
as comparison (solid line) with the median vertical error bar in the left corner. A 
vertical arrow signs the reference galaxy used for the theoretical calibration. 



9 r r 




2 3 4 5 

log (pc) 

Fig. 6. As in Fig.® with d = 0.6 calibrated with M%/L° = 15, 10, 5, 1 (from top to 
down, with the same previous line types) and the fit of Djorgovski & Davies (1987) 
(solid line), in Lick re band, as comparison with the median vertical error bar in 
the left corner. A vertical arrow signs the reference galaxy used for the theoretical 
calibration. 

10; Xe = (1.70)2; Fi(Xi) = 2), we obtain ki = 3.79 ■ 10^ Lq / {KpcY ; = 
1.31 ■ LQ/{Kpcf and the total luminosity : L° = 3.78 ■ lO^L©. 
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To transform the photometric data of our reference galaxy from the JohnsoiT^ 
BVR system into Lick re band pass we use the Djogovski (1985) transforma- 
tions: (Vj — tg) = 0.031 + 0.681(5 — V),j, assuming a mean {B — V)j ~ 1, 
(^g)© = 4.43 and log/g = — 0.4((/i)e) — 26). At this level of first approximation 
comparison we do not consider how could change re and cTq of our reference 
galaxy by changing the photometric color band. 

The comparison looks fairly well. Inside the vertical median bar the edge-on 
observed FP of DD87 lies between the two theoretical FPs corresponding to 
M^/L° = 10 and M^/L° = 5, in the case d = 0.5 (FigED. The theoretical 
values A*, B*, —OAB* /A* are 1.33, —0.83, 0.25 and the observed ones are: A = 
1.39 ± 0.15, B = -0.90 ± 0.10, -OAB/A = 0.26 ± 0.06. Then inside the error 
bars the two results coincide even if the theoretical straight lines turn out to be 
a little bit steeper in respect to the observed one. The other case d = 0.6 (Figl6]) 
looks better from the slopes point of view: the theoretical values become: 
A* = 1.43,5* = -0.86, -0.45* /A* = 0.24 in very good agreement with the 
fit result. The two slopes, theoretical and observed, are then nearer, even if an 
higher value of M^/L° = 15 seems to be preferred without exclusion of the 
case M^/L° = 10. Concluding, from this first approximation comparison both 
cases d = 0.5, d = 0.6,M^/L° = 10 are acceptable inside the error bars of 
observed fit with an M'^/L° = 5 to be preferred in the case d = 0.5 from the 
reference galaxy (vertical arrow) forwards. On the contrary the case d = 0.6 
seems to request an higher ratio M'^/L° at least of 10. 



8 The most relevant way 



The most relevant way to understand the physical meaning of FP occurs as 
soon as we wish to translate the theoretical FP described by the relationship 
fl62l) into the K-space (Bender et al. 1992). Now we know the second relation 
fl63l) as an equation by King's model, that is Eq. fl59l) . The fl63l) tells us how Jg 
has to scale with M/? and Mb- We know already it from LSI: 



~ Mi,; i = I = 2 ^' f (73) 



y(3-d) 



but it is Eq. fl59l) which allows us to understand deeply why this scaling law 
has to be followed. We can indeed to recover the fl73|) by Eq. fl59|) as soon as 



^ Actually we used the UBVRI bands for the combined Johnson-Cousins-Glass 
system and the solar corresponding values as given by Binney &: Merrifield (1998, 
Chap. 2, pg.53). 
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-(^-1) 



we remember that: 

PoD-^M„"'' " (74) 
according to subset. 6.1. By using Eq. fHSj) it follows: 



yr'^MO'D (75) 



1 



where from cosmology (Eq. HOl) . a/) ~ M^' . 
Remembering that: 



l-d 



kj^lkM = ^ = Ml- = M^'-^ (76) 
Mb 

we obtain, at fixed Cb'- 



h ~ r^PoD {—) ' 1^ ~ ■ M" ^-'^ ■ Mi' (77) 



-(4--l| d_ - <* l-d 



with the exponents 



( ?, \ d d d I I 
x = -[--l]+- + — — - - - — — + - 78 

\j ) J 6 — d J 3 — d y 

o y-(3-^) ^ 

Y (3 - d) 

1 C? 1 — d ^ ^^g^ 

3 — (i 3 — c? 3 — 
which matchs exactly (173ll . 

The physical meaning is the following: light does not follow the visible matter 
because it depends, via stars, on the deepness of the gravitational potential 
well which is determined by the two central densities of DM and baryons linked 
together by the equipartion of virial energy, that is by the Eq. flST]) . 



8.1 If.- scaling 



If we take into account two galaxies corresponding to two different baryonic 
masses (M^,Mb) but characterized by the same m = M^/Mb = M^/M^, 
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Cb and Cd (i.e., the same CVM, xm = XM{d,m,CB,CD)), h has to scale as: 



= {Mb/M^^Y (80) 

which reduces to: 

IJi: = {Mb/M"^)-^ (81) 
if we are in the typical mass range of M^) ~ IO^^Mq, that is 7' ^ 2. 

8.2 Te yCTo- scaling 

From Eq. flB^ we know that: 

.2 - ^ (82) 

reC2C3 

where Ve scales in turn with Mb as follows: 

r, ~ M'J'k-J^' (83) 

due to the King's model relationship, /cA/r^ ~ Mb, at fixed Cb- But at the 
Clausius' virial maximum configuration, has to scale as (LSI): 

r, ~ m''M|; r = J/2; R = I/7' (84) 

For the two galaxies considered before, at fixed xm, it follows that: 

rjr: = {MB/M°j,f^' (85) 
as soon as ~ M^^'^^^^"' , and: 

{ao/Kf = (Mb/M^s)'^ (86) 
9 Theoretical tilt equation in K-space 



Following Bender et al.(1992) we have to build up the theoretical tilt equation 
in the K-space. From Eqs.f l80|l85]|86|) we obtain for ni and k,^, respectively: 



Ki = (log (To + logr, 
= (log«)2 + logr:) + log(MB/M°))/v^ 
log(MB/M°) = V2k, - [log«)2 + logr^] 



^7) 



23 



and 



^3 = {logal - log/e - logre)/v^ (89) 
= (log«)2 - log/: - logr:) + ■ log(MB/M^))/v^ 

7' - 2 

ae = Qeirirec, d) = ( 1(7', d)) (90) 

7 

Inserting the link (Eq.( [88])) between (ki, K3), we obtain the tilt equation: 



.3 = ■ + - - ^+ - (91) 

The last equation tells us how the FP becomes degenerate in respect to the 
cosmology. Indeed, as soon as we are in the typical galaxy mass range, so that, 
y ~ 2 (Gunn, 1987, Silk, 1999, Chap.3), then -/ ^ a, according to Eq.(l8T1). 
and ae —>■ a which is depending only on d, i.e., on the halo mass distribution. 
In this range the galaxy FP becomes: 

fTj:; , (1 - a) log«)^ -jl + a) logr° - log J° 1 - ^ .qo^ 

K3 = KiH -j= ; a = j—^ (92) 



The last equation shows that: ifd=l^a = and then ^3 = const., which 
means the tilt disappears! 



9. 1 Calibration 



Again as in subsect. 7.3, we calibrate the FP by the elliptical galaxy PGC045032 
(1300.4 + 2807) of Coma cluster {zc = 0.023) which has been fitted by King's 
profile (Oemler, 1976), characterized by: 

Cb = 10; Fizo) = 0.7294; = (1.70)^; Fi(Xt) = 2 (93) 

We choose for all the galaxies of the theoretical plane the same Clausius' virial 
maximum configuration, corresponding to: 



Cn = CB = 10; m = 10; d = 0.5; Uy = 0.0697; uns = 0.9039 (94) 
(iyD)M = 0.5463; Mm = 0.0103; xm = 0.9; 

If our reference galaxy has a ratio: k^/kM = 1/5, B — V ^ 1, we obtain, in B 
band: 
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log«)^ = 3.83 (95) 
\ogI° = 1.98(-^ 1° = 95.2Lq/pc'^; L° = 2.61 ■ W^Lq) 

logr° = 0.33 (-^ r° = 2.UKpc) 

= 0.88, k1 = 2.94 (96) 

Then the theoretical ti/t equation (l92l) in k- space ( d = 0.5 ^ a = 0.2; 7' ^ 2) 
becomes: 

= O.I6K1 + 0.40 (97) 
to be compared with that given in B band by Burstein et ah (1997): 

4 = 0.15^ + 0.36 (98) 

The two equations are plotted in FigjTl The difference in K3 at fixed ki (Figl?]) 
turns out to be: An^ = 0.07, a little bit over the assigned FP tightness: 
a{K3) = 0.05. 



1.5 

S 1 
0.5 



Kj [Mass] 

Fig. 7. Comparison between the theoretical tilt-equation (Eq. (j97p : dashed line) 
(d=0.5) and that derived by the observational fit in B-band of Burstein et al.(1997) 
(solid line). The discrepancy between the two straight lines is AK3 = 0.07 at the 
reference galaxy used for calibration (signed by a filled square) less than 2 times 
the assigned tightness of FP, a{K3) = 0.05, shown on the left corner. 

It should be noted that the DM halo of our reference galaxy turns out to be 
characterized by: 

Md = mM°B = 1.3 ■ lO^Mo; 7' ^ 2; (99) 

Its formation redshift has been calculated better than by Eq. ( H2i) . using the 
subroutine of Navarro et al. (1997) in which zp is precisely defined as the z at 
which half of final mass is in progenitors more massive than 1% of the final 
mass. At this zp = 3.2 {h = 0.75) the corresponding overdensity (in units of 
the critical density at ^ = 0) becomes Sc = 2.42 ■ 10^ which corresponds to 
the value of central density: {pod)nfw — 3.3 ■ 10~^Mq/pc^, for NFW-profile 
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Fig. 8. Comparison between the theoretical tilt-equation (Eq. (j92p : dashed line) 
(d = 0.6; K3 = 0.14ki + 0.46) and that derived by the observational fit in B-band of 
Burstein et al.(1997) (solid line). The discrepancy between the two straight lines is 
Ak^ = 0.06 at the reference galaxy used for calibration (signed by a filled triangle, 
Ki = 2.91; Kg = 0.85), about the assigned tightness of FP, (t(k3) = 0.05 shown on 
the left corner. The limit case d = 1 for which the tilt desappears, is also plotted 
(dotted line) with the relative position of the calibration galaxy (filled circle). 

used by Navarro et al. (1997). Rescaling this value to the cored power-law 
profile of Eg. (1291) {d = 0.5) by a factor of about 1.17, to match the theoretical 
value of /e given by Eq. fl59|) we need of contraction factor c ~ 1/15. That 
appears consistent with the limits found by Klar & Muecket (2008) (sect. 6), 
considering that in this context most of baryonic matter is in stars. 

9.2 Discussion and conclusion 

Both the comparisons of theoretical FPs either with that of Diorgovski & 
Davies (1987) in Lick or that of Burstein et al.(1997) in B band result very 
satisfactory. Moreover the edge-on theoretical representation, [k3,ki), of FP 
(FiglH]) appears to rotate around the calibration point as soon as d increases 
according to Eq.( !92l) . It becomes horizontal when d becomes equal 1. In this 
limiting case: d = 1 — a = 0; = Kg whichever is ki and then the FP looses 
the tilt (FiglH]) as already underlined in the previous papers LSI, LS5. 

The new order of approximation of TCV by the inclusion of a King-like B 
component allows us to understand more deeply the physical reason of the FP 
tilt. The CV theory is born to try to explain the FP (firstly of ETGs) without 
breaking the homology +virial equilibrium. When the galaxy system is looked 

^ For this profile, central density means density at about one half of scale radius. 
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as a whole system a dynamical explanation of the observed tilt turns out to be 
impossible without breaking homology: the exponents A, B become rigorously 
2,-1 due to the virial equilibrium. The only way to change the exponents 
is to split the system into two subsystems: one of baryons, the other one of 
DM as done by the tensor virial theorem. Then the double system may gain 
a new symmetry due to the equipartition of Clausius' virial energy between 
dark and visible matter. When that occours, the virial configuration becomes 
special because its CV energy reaches a maximum value, i.e., a minimum of 
its random kinetic energy to obtain equilibrium. This allows to the exponents 
of cTq and Je to decrease their absolute value if the bulk of baryonic matter lies 
inside a dark matter halo distribution like a power law pn ~ r"'^; d ^ 0.5. But 
that is not enough. The reason of the tilt, at this higher order of approximation 
of TCV, appears more clearly. It lies in that: light does not follow the visible 
matter because the star formation is due to the two potential wells, of baryons 
and of DM. Their depths are linked together by the equipartition of virial 
energy, that is the relationship given by Eq. flSTl) . Then the main conclusion is 
the following: in our approach the galaxy tilt is neither due to different DM 
fraction which enters into the dynamical mass of Eq. fHOl) . in order to increase 
the observed ratio M^yn/ L (i.e., K3) at increasing M^yn (i.e., ki) starting from 
a fixed mass-luminosity ratio T^, = Mb/L for all galaxies. That has been 
already understood by Ciotti et al. , since 1996, who realize a fine-tuning 
was invoked. Nor it may be explained by trying to tune Mdyn/L by different 
DM amount, again with a constant T.,,, assuming the galaxies are located on 
the CV maximum configuration as Valentinuzzi (2006) tried whithout success 
(even if his approach was substantially different in respect to the TCV). In 
our approach the fraction of Mdyn in DM is always constant, as soon as we 
locate the ETC at a fixed amximum: Xm = XM{d, m,CD,CB)- That is: 



In this way the: 



That means the DM amount which enters into M^yn increases at increasing Mb 
but in a way proportional to Mb, so that any tilt is produced if we start with 
T* = const.. The TCV offers the dynamical mechanism in order to change T,,, 
exactly ~ Mg. 

The other important answer of TCV is why the FP as a whole appears to be 
degenerate in respect to cosmological density perturbation spectrum. Looking 
at Eq. (!9T]) the tilt looks apparently not independent of cosmology because 
is depending explicitely by 7' which appears also inside the exponent I. But in 
CDM scenario the effective index on typical galaxy scale, for a scale-invariant 
initial spectrum, is indeed approximately —2 (Silk, 1999, Chap. 3) and that 




(100) 



27 



corresponds to 7' ~ 2. So the Ue degenerates into a which is independent of 
cosmology. It depends only on the DM distribution d. That is also the reason 
why the ratio L/Mb is totally independent of the cosmic perturbation spec- 
trum and of mass ratio m as we have already proved in LSI. The degeneracy 
is broken as soon as we look at the projections into the coordinate planes, as 
shown in LSI. 

Many problems are still open. We have considered a unique maximum for all 
ETG and a unique Cb for the King-like B component. What occurs by chang- 
ing them. What about the other parameters involved? What happens as soon 
as we change 7'? Moreover we may wonder how the results do change moving 
to a ACDM scenario, even if we expect no significant variations essentially 
because the mass variance does not change too much in this last cosmology. 
But the main point is: why has the DM distribution which contains the bulk 
of baryons follow a density power law of the kind pu ~ r~'^;d = 0.5 -r 0.6 
in order to produce the observed tilt. Many efforts have been devoted to this 
problem. Theoretical arguments based on dynamics (Miicket & Hoeft, 2003) 
and thermodynamics (Secco et al., 2007) strengthened by observations lead 
on this direction even if numerical simulations seem to prefer ~ 1 (see, e.g., 
Bindoni, 2008). At the moment a definitive answer to this crucial point does 
not exist. 
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A Appendix 

A.l King's dimensionless density profile 

The King (1962) spatial density profile of Ea U6l in the explicit form, is: 



vrrc [l + in/rcf] 

\ h + jrt/rcf 
■ I [ 1 + {r/r,f 

If we define = r/ciB = r/rt and Cb = rt/rc, which means Cb£,b = ''/''ci we can write: 



1 + {rt/rcY 
1 + (r/rcf 



1 + {r/raf 
1 + {n/rcY 



(A.l) 



{rt/r^Y - {r/vaY 
l + (rt/re)2 



28 



ktvi 



1 + {CB^sr 



l + {Ci3^i3)2^ ^ 



1 + ^1 

1 + (Cs5s)2 



1 + C| 



1 + C| 



(A.2) 



and: 



Po = p(« = l/Cs) 



l + Cl 



1 + Cl 



i + cl 



i + cl 



(A.3) 



Then, the King's profile, normalized to the scale radius density value, is: 



/B(eB) 



l + (Cses)2 



l + (CBiBf 

1 + C% 



Po 1 + {CsCs)2 

cl-(CBiBr' 



1 + C| 



(A.4) 



where: 



H ■- 
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1 
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(A.5) 



A.2 DM dimensionless density profile 



The cored power law that describes the DM density profile is: 



PD{r) ■ 



2poD 



1 + {tItodY 



(A.6) 



where ToD is the scale radius and poo the density value at the scale radius. In the usual way, once defined 
= r/ajj and Co = ao/roD, which means Cu£,d = ^hoD-, we can write: 



Pn{(,n) ■■ 



2poD 



1 + (Cd^dY 
Then the normalized DM profile becomes; 



(A.7) 



PdUd) ^ 2 

PoD 1 + {Cd^dY 



(A.8) 
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A. 3 Calculation of Clausius Virial 



We define all the coefficients we need for the Clausius trace. All of these are depending on different values 
of Cs, Ci3, d, m = Mjy/Mg and the computation of them, where it was not possible in an analithical way, 
was performed numerically by the software Derive for Windows 6.0. In the Tabl lA.llA.2l IA.3I I are listed the 
values of all the coefficents for different parameters. Here we present how they were calulated, for example 
in the case d = 0.5: 



/o 

»1 



H 



l + iCB^BY 



1 



1 + {Cb^b)^ 
1 + Cl 



{Cb^b? 



Fb{^b) = 2 / fB{^B)iiBdiB 



1 + C| 



l + {Ci3Ci3)2 



1 + (Cb5b)^ 
1 + [Cb£.b? '^^ 



H 



[z{iB)]l 



cl - [Cb^b? 



Cl 



■ ^Bd^B 



(A. 10) 



where: 



ZiiB)=--^ 



Cl + 1 



cl V (CBCfl)2 + l 



l + (CsCs)2 



1 



cl 



^In [(Cfl5fl)2 + l] + 



^2 
^B 



In 



CBy/Cl + l-^l-^B + Cl + l 



(A.ll) 



I^QB 



16 {1^3)1, 



Fl{iB)diB 



(A.12) 



FD{iD) = 2 I fD{^D)iDd£,D 

ho 

"1 



xO.5 



e„ 1 + (CD£.Df 

iHy'Coe.D + l) , 4(Cpgp + 3)^(7ogg 2go 

Cd 



Cl 



3Cl 



in 



8 ln( y^CdCd + 1) 81n(vC^J+l) 



Cl 



2(CDiD + 3) - VC^{3^/Cd^d + 2Cd - S^C^ + 6) 



3Cl 



(A.13) 
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Jo 



15C3 



Co C2 



(A.14) 



12C|, - IbC^J'^ + 20Cd - 30C]^'^ + 60 12 In + 1 



5C' 



5/2 



'^D 1 + ^Co^n 



(A.15) 



If we define x = aB/ao we can exprime (jj in terms of in the following way: = ^b (f^) = (b^- 



Jo 



Wext{x)= I FB{iD)^^^^-iDdiD 



din 



(A.16) 



i/v ■ 



8 (i's)m (''d)m 



• m ■ Wext (x) 



(A.17) 



At the end, Clausius Virial, normalized by the factor GMBF/ao can be expressed by: 



Vb (x) = - 



(A.18) 
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Table A.l 

Physical parameters of King's models as function of concentration, Cb- 
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Table A.2 

Physical parameters of the DM halo described by a cored power law density profile 
with exponent d. 
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